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ABSTRACT 

The microquasar SS 433 is an interacting massive binary consisting of an 
evolved mass donor and a compact companion that ejects relativistic jets. The 
mass donor was previously identified through spectroscopic observations of ab- 
sorption lines in the blue part of the spectrum that showed Doppler shifts associ- 
ated with orbital motion and strength variations related to the orbital modulation 
of the star-to-disk fiux ratio and to disk obscuration. However, subsequent ob- 
servations revealed other absorption features that lacked these properties and 
that were probably formed in the disk gas outfiow. We present here follow-up 
observations of SS 433 at orbital and precession phases identical to those from 
several previous studies with the goals of confirming the detection of the mass 
donor spectrum and providing more reliable masses for the two system compo- 
nents. We show that the absorption features present as well as those previously 
observed almost certainly belong to the mass donor star, and we find revised 
masses of 12.3 ± 3.3 and 4.3 ± O.SMq for the mass donor and compact object, 
respectively. 
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INTRODUCTION 



The uniqu e system SS 433 is an X-ray binary star that falls in the class of microquasars 



( iFabrikal 12004 ) . The system exhibits relativistic jets originating from a compact object 
surrounded by an accretion disk. The companion star in the system contributes only a 
small fraction of the total light and is thus difficult to detect. Because the accretion disk 
and its wind are the dominant source of light in the system, reasonabl e measurements have 
been made of the radial velocity amplitude of the compact object (e.g. iFabrika fc Bychkova 
1990l : iGies et al.l l2002al ). Models of th e precessing jets and eclips es have also provided a 
very accurate system inclination (78?8. iMargon fc Anderson! Il989l ). Therefore, identifying 
the mass donor and accurately determining its radial velocity amplitude will lead to direct 
kinematical masses for both components. 



Th e first potential observations of the mass donor st ar were reported bvlGies. Huang, fc McSwain 
(2002b) with additional observations published later by Cherepashchuk et al. ( 20051 ). These 
studies detected faint absorption lines present in blue spectra of SS 433 that exhibited the 
Doppler shifts expected for the donor star. Finding these lines is difficult, not only due to 
the small relative flux contribution of the A-star, but also to the active emission spectrum 
from the accretion disk, jets, and strong disk wind in the system. The most convincin g 
observations to date of the mass donor spectrum were obtained by iHillwig et al.l (120041 ). 
Their spectra suggested that the donor star has a spectral classiflcation of A3-7 I, and 
their radial velocity measurements led to component masses of Mq = 10.9 ±3.1 Mq and 
Mx = 2.9 ± 0.7 Mq (where the subscripts represent the optical companion star and the 
X-ray emitting compact object X). 



This identiflca. tion o f the mass donor star w a s cha llenged by I Charles et al.l (120041 ). 
Barnes et al.l (120061 ) . and IClark. Barnes, fc Charles! (12007! ). These authors present spectra 
of SS 433 that show features found in an A-supergiant, but that do not follow a coherent 
orbital radial velocity curve. Their observations show a larger scatter, were performed over 
a greater variety of orbit al and precess i onal p h ases, and have a mu ch lo wer average yelocity 



than those presented bv IHillwig et al.! (12004! ). iBarnes et al.l (120061 ) and IClark et al.! (12007! ) 
suggest that these A-supergiant absorption lines arise in an accretion-driven outflow. Thus, 
their work demonstrates that not all the absorption features in the spectrum of SS 433 form 
in the photosphere of the mass donor star. 

Some of the differences between these investigations are due to the timing of the ob- 
servations and the inherent time variability of the spectrum. Given the rela tive faintness 



of the mass donor star compared to the super-Eddington disk in the system, IHillwig et al. 



( I2OO4J ) planned observations for orbital phases when the disk is partially eclipsed and for 
precessional phases when the donor star appears above the disk plane at orbital inferior 
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conjunction. This strategy maximizes the opportunity to reduc e the glare of t he di sk and 
avoid obscuration by the extensive disk gas. On the other hand, iBarnes et aLl (120061 ) made 
observations at a variety of other orbital and precessional phases and made plausibility ar- 
guments that the donor star spectrum might appear at some of these (but was not in fact 
detected) . 



Given the need to r e solve the differing interpretations in the papers of iHillwig et al. 



(120041 ) and 



Barnes et a. 



(120061) . we undertook follow-up observations to repeat the obser- 



vations of iHillwig et al.l (j2004j ). If the absorption features identified do originate in a disk 
outflow, then the highly variable nature of the outflow in SS 433 suggests that repeat obser- 
vations will show significant scatter from prior observations, if in fact the absorption features 
are still present. However, if the absorption lines originate from the mass donor star, then 
the highly stabl e and periodic natu re of the orbit should produce observations consistent 
with the data of IHillwig et~aD tooj ). 



Observations and Reductions 



Observations were acquired with CMOS on the Gemini North Telescope on UT 2006 
June 7-13. The resulting spectra cover the wavelength range 3767-5227 A with a dispersion 
of 0.235 A pixel"^ and resolution R{X/AX) = 9580. Observations were made for approxi- 
mately four hours each night, and all the individual spectra from each night were coadded to 
improve the S/N. The coadded spectra were shifted to a heliocentric frame and continuum 
normalized by fitting regions free of strong emission lines. 

The observations were t imed to corre s pond to mid-eclipse and greatest disk opening 
angle, as with the data from IHillwig et al.l (120041 ). The spectra cover orbital phase 0.85 < 
4) < 0.32 and precess ional phase 0.02 < if) < 0.06, according to the orbital ephemeris of 
Goranskii et"al] Jl998l ). 

HJD 2, 450, 023.62 + 13.08211E 



and, for disk precession, the model ephemeris of iGies et al.l (l2002al ) 



HJD 2,451,458.12 



162.15E. 



The orbital phase coverage overlaps with and extends that of 



Hillwig et al, 



mm . 
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3. Spectral Appearance and Radial Velocity Analysis 

The resulting spectra from seven consecutive nights around mid-eclipse show the ex- 
pected broad "stationary" emission lines. The brightest emission lines of hydrogen and 
helium all showed strong P Cygni absorption components every night except on the fourth 
and last nights of the run. The B.j3~ jet line was present in the spectra, but weak. 

The new series of spectra appear to show considerably more emiss ion contamination i n 



the range AA4500-4625 A than seen at the time of the investigation by lHillwig et al.l ( 20041 ). 
and consequently, it was not possible to measure unambiguously the same set of absorption 
lines they analyzed. The only clearly problem free regions larger than ~ 25 A in the observed 
spectral range are the intervals of AA4750-4830 A and AA4950-4985 A. These regions show 
clear absorption features in the spectra. An interesting case is that of night seven, in which 
the only visible absorption features are interstellar in origin. Their clear appearance in that 
spectrum allowed us to subtract the interstellar absorption from the spectra obtained on the 
first six nights. This was especially helpful near 4762 and 4780 A where interstellar lines 
overlapped apparent stellar absorption features. 

The resulting, interstellar-subtracted, and normalized spectra of SS 433 in the wave- 
length region described above are shown in Figure [T] along with a spectrum of an A- 
supergiant, HD 9233 (A4 lab). The HD 9233 spectrum was obtained in 2005 August with the 
medium resolution spectrograph (MRS) on the University of Texas Hobby-Eberly Telescope 
(HET). It has been smoothed to match the resolution of the SS 433 spectra and intensity 
scaled to match the mid-eclipse line depth as described below. A comparison of the inter- 
stellar spectrum of SS 433 from night seven with that of HD 9233 shows that the latter has 
no appreciable interstellar absorption. We confirm that the the lines of an A-type supergiant 
absorption spectrum are once again present in the spectrum of SS 433. 

Line identifications for major absorption lines are also shown in Figure [1] All the 
spectra in the figure are Doppler-shifted to the heliocentric rest frame. The SS 433 spectra 
were shifted using the radial velocity values discussed below. The HD 9233 spectrum was 
shifted using a measured heliocentric radial velocity of —69.5 ± 1 . km s~^, significantly 



different than the value of —34 ± 2 km s ^ given by iHillwig et al.l (120041 ). It is likely that 
HD 9233 is a radial velocity variable. 

We measured radial velocities for the ab sorption lines in the SS 433 spectra using the 



same method applied by lHillwig et al.l (120041 ). i.e., by calculating cross-correlating functions 
(CCFs) of the absorption line regions in each SS 433 spectrum with those in the spectrum of 
HD 9233. The spectrum of HD 9233 was first rescaled in intensity to match t he absorption 



features of the mid-eclipse spectrum of SS 433. The scaling factor found in iHillwig et al. 



- 5 - 




Fig. 1. — SS 433 spectra obtained with GMOS at Gemini North in 2006 June along with 
a spectrum of the A4 lab star HD 9233 {bottom) obtained with the MRS on HET in 2005 
August. The HD 9233 spectrum has been smoothed and scaled to match the SS 433 spectra 
(see text). Spectra have been offset for clarity. The corresponding orbital phase for each 
spectrum is shown on the right side of the plot. 
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( 120041 ) of 0.36 ± 0.07 was also found to be appropriate for this data. We then added the 
radial velocity of HD 9233 to the relative velocities from cross-correlation to transform them 
to an absolute scale. The resulting heliocentric radial velocities for SS 433 are given in Table 
[Hand are plotted as a function of orbita l phase in Figure [2] (together with the radial velocity 
measurements from lHillwig et al.ll2004l for comparison). 



The CCF fitting was adversely affected by the emission near the Cr II A4824 line so that 
line was removed from the CCF fitting region. Also, the fit on night six is clearly dominated 
by the strongest line (Fe I A4957). To determine if this line might be consistently offset from 
the full sample CCF values, and thus offset the night six radial velocity value, a CCF fit 
was run using just this line. The resulting radial velocity and CCF amplitude results for 
nights 1-6 were well within the individual error ranges of the full CCF fits and there were 
no apparent systematic differences. 

The absorption lines shown in Figure [1] appear to be strongest near mid-echpse (phase 
0.0) when the flux dilution from the disk light is minimized. This variation in strength is 
represented in Figure B| by the relative CCF amplitudes for the 2006 June Gemini North data, 
along with those from lHillwig et al.l (120041 ). The solid line in this figure shows the predicted 
variation in strength for a constant flux component that is diluted by the phase- variable flux 
co ntribution from t he disk. Th is function was derived from the sch ematic i?-band light curve 
of lGoranskii et al.l (Il997l . ll998l ) as described in lHillwig et al.l (120041 ). Since the disk flux varies 
stochastically, the actual eclipse light curve during the time of the 2006 June observations 
may have differed from this prediction. Nevertheless, it appears that the absorption lines 
did strengthen and weaken as expected during this particular eclipse. 



Table 1. Absorption Line Radial Velocities 



Date 






K 




(HJD - 2,450,000) 






(km s~^) 


(km s-i) 


3893.99 


0.022 


0.852 


25 


4 


3895.04 


0.029 


0.932 


57 


4 


3895.96 


0.034 


0.003 


77 


4 


3897.02 


0.041 


0.084 


96 


4 


3898.01 


0.047 


0.159 


125 


6 


3899.03 


0.053 


0.237 


135 


6 


3900.03 


0.060 


0.314 
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Fig. 2. — The phase-folded radial velocity curve for the absorption spectrum in SS 433. 
Shown are t h e Ge mini North data from 2006 June [filled circles) and the data from 



Hillwig et al.l (120041 ) (hollow squares). The dashed lines show the ±2a semiamphtude so- 



lutions relative to the calculated best fit. 
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Fig. 3. — The normalized CCF amplitude vs. orbital phase curve for the SS 433 absorption 
spectrum. Shown a re the results from th e Gemini North data from 2006 June {filled circles) , 
the CCF data from lHillwig et al.l (120041 ) {hollow squares), and the prediction {solid line) for 
a constant flux photospheric spectrum diluted by the variable flux from the accretion disk 
(according to the mean i?-band light curve). 
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Discussion 



Hillwig et al.l (120041 ) discuss three criteria that must be met if the absorption line pattern 
is to be rehab ly associated with the supergiant: (1) The radial velocity variation with orbital 
phase should appear with the ascending branch centered at orbital phase 0.0 (i.e., at inferior 
conjunction at the time of the optical and X-ray mid-eclipse). Furthermore, the systemic 
velocity should be approximately the same as that of the nebula (+56 to +75 km s^^; 
Boumis et al.ll2007l : iLockman. Blundell. fc Gossll2007l ). The new radial velocity data (Fig. [2]) 
agree with these predictions and with the earlier radial velocity measurements. (2) The 
absorption line depths will usually be weak and modulated in strength by the orbital curve 
and precessional phase. The observed spectrum is a composite of disk and stellar light, so 
that the stellar lines will always be diluted by the continuum flux from the disk. Furthermore, 
the accretion disk is probably optically thick and vertically extended, so the supergiant will 
often be hidden when the star is behind the disk plane from our line of sight. The new 
sp ectra show th e same variation in absorption line strength (Fig. [3]) as found previously 
by iHillwig et al.l (120041 ). consistent with our expectations for the varying supergiant to disk 
flux ratio. (3) The presence of a super-Eddington disk in SS 433 implies that the system is 
undergoing rapid mass transfer from an evolved and Roche-filling mass donor. The fact that 
the absorption line spectrum resembles that of an A-supergiant with a projected rotational 
l ine broadening tha t matches expectations for a Roche-filling star in synchronous rotation 
( IHillwig et al.ll2004l ) is fully consistent with the mass transfer picture. 



All three criteria are met in the new ob servations, and the r emarkably consistent agree- 
ment between the new results and those of IHillwig et al.l (120041 ) strongly support the con- 
clusion that the absorption pattern in the spectrum of SS 433 observed near mid-eclipse and 
precessional phase zero is that of the mass donor star. On the other hand, the behavior of an 
absorption spectral component associated with a disk or disk wind will be highly tirne vari - 
able , and we suspect t hat the kind of absorption features described by iBarnes et al.l (120061 ) 
and I Clark et al.l (120071 ) do originate in a disk outflow as they propose. 



The radial velocity curve also provides us with the opportunity to revise the kine- 
matical mass calculations for each component. The best-fit radial velocity curve for the 
combined sets of data shown in Figure [2] gives a semiamplitude for the mass donor star of 
Kq = 58.2 ± 3.1 km s~^ and a systemic velocity of 7 = 73 + 2 km s~^. A dopting a compact 
object semiam plitude of Kx = 16 8 + 18 km s~^ 



see 



Hillwig et al.ll2004j ) and system incli- 
nation of 78?8 (IMargon &: Andersonlll989l ) gives component masses of Mq = 12.3 + 3.3 Mq 
and Mx = 4.3 + 0.8 Mq. The larger revised mass for the compact object establishes it 
much more firmly as a black hole candid ate. The derived mass ratio Mx i Mo = 0.35 



is consistent with limits from the optical (lAntokhina &: CherepashchukI Il987l ) and X-ray 
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( lAntokhina. Seifina. fc Cherepashchuklll992l : iFilippova et al.ll2006l ) light curves. 

We caution that our estimate of Kq may differ from the actual geometric value if 
the inner hemisphere of th e A-supergiant is significantly heated by flux from the disk 
( Cherepashchuk et al.l [20051 ). The amount of irradiation present is very uncertain since the 
disk vertical extensions may be sufficient to block most of the high energy radiation from 
the central accretion zone from reaching the facing hemisphere of the supergiant. Further- 
more, the effects of such heating on line formation are complicated. For example, a hotter, 
brighter inner hemisphere will have a center of light that is shifted towards the companion 
(so that the measured radial velocity shifts would be less than actual). On the other hand, 
if the temperature is high enough in the irradiated hemisphere to shift the plasma to higher 
ionization stages or to reduce the temperature gradient in the atmosphere, then the lines we 
observe would be weaker or absent in the surface elements of the inner hemisphere (shifting 
the measured radial velocities to larger than actual). The clues from the available spectra 
suggest that irradiation probably plays a minor role since the decline in CCF strength fol- 
lowing mid-echpse follows approximately the expected curve for flux dilution alone (Fig. [3]). 
However, the absorption lines do weaken and disappear in the final spectrum obtained near 
orbital quadrature phase, and this behavior may be due to some contribution of irradiation 
or the disappearance of the supergiant behind an optically thick accretion disk. If irradiation 
of the inner hemisphere were to change the ionization stage, however, then at quadrature 
the absorption lines would be visible from only half of the projected stellar surface, meaning 
the line strengths should be no greater than half that predicted by the eclipse model {solid 
line) in Figure O 



5. Conclusions 

We have shown that the behavior of the A-supergiant absorption spectrum observed in 
SS 433 during eclipse and near precessional phase zero is repeatable. This suggests that the 
origin of the spectrum is a consistent, well-behaved source. We suggest that the source is 
the mass donor star in the system rather than the disk or a disk wind. It is interesting that 
SS 433 would show such similar spectra from both the mass donor star and the disk/disk 
wind. However, the appearance of such absorption might be anticipated if the disk outflow 
material is clumpy and has a temperature in the line forming region comparable to that of 
the supergiant. 

The radial velocity curve resulting from observations of the mass donor star then allow 
calculations of the kinematical masses of the two components. We find masses of Mq = 
12.3 ±3.3 Mq and Mx = 4.3 ±0.8 Mq for the mass donor and compact object. It is possible 
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that the inner hemisphere of the mass donor is irradiated, which would change the mass of 
both components, potentially reducing the compact object mass to below the neutron star 
limit. The behavior of the CCF amplitude suggests that a strong irradiation effect is not 
present in the data presented here. If this is the case, then the compact object in SS 433 
can be firmly identified as a black hole candidate. 

We thank the staff of the Gemini North Observatory for their help in obtaining the 
GMOS spectra described here. This material is based upon work supported by the Na- 
tional Science Foundation under Grant No. AST-0606861 (TCH). Institutional support has 
been provided from the GSU College of Arts and Sciences and from the Research Program 
Enhancement fund of the Board of Regents of the University System of Georgia, admin- 
istered through the GSU Office of the Vice President for Research (DRG). We gratefully 
acknowledge all of this support. 
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